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Mid-rapidity open charm spectra from direct reconstruction of D°(D°)~^ K t tt in d+Au colli- 
sions and indirect electron/positron measurements via charm semileptonic decays in p+p and d+Au 
collisions at ^/s NN = 200 GeV are reported. The D°(D°) spectrum covers a transverse momentum 
(pr) range of 0.1 < pr < 3 GeV/c whereas the electron spectra cover a range of 1 < pt < 4 
GeV/c. The electron spectra show approximate binary collision scaling between p+p and d+Au 
collisions. From these two independent analyses, the differential cross section per nucleon-nucleon 
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binary interaction at mid-rapidity for open charm production from d+Au collisions at RHIC is 
do" ( ^ iv /dj/=0.30±0.04 (stat.)±0.09(syst.) mb. The results are compared to theoretical calculations. 
Implications for charmoniumm results in A+A collisions are discussed. 

PACS numbers: 25.75.Dw, 13.20.Fc, 13.25.Ft, 24.85.+p 



Hadxons with heavy flavor are unique tools for 
studying the strong interaction described by Quan- 
tum Chromodynamics (QCD). Due to the large mass 
of the charm quark (~ 1.5 GeV/c 2 ), charm quark 
production can be evaluated by perturbative QCD 
(pQCD) even at low momentum through the intro- 
duction of additional scales related to the charm 
quark mass 0, 0- Therefore, theoretical calcula- 
tion of charm hadron total cross section integrated 
over momentum space is expected to be less affected 
by non-perturbative soft processes and hadroniza- 
tion Q . Systematic studies of charm production in 
p+p and p+nucleus collisions have been proposed 
as a sensitive way to measure the parton distribu- 
tion function in nucleons, and nuclear shadowing ef- 
fects E|. At RHIC energies, heavy quark energy 
loss [J^ charm quark coalescence ja IE S, possi- 
ble J/ij) suppression 01 i an d charm flow |ll| have 
been proposed as important tools in studying the 
properties of matter created in heavy ion collisions. 

Identification of charmed hadrons is difficult due 
to their short lifetime (ct(D°) = 124 ^m), low pro- 
duction rates, and large combinatorial background. 
Most measurements of the total charm cross section 
in hadron-hadron collisions have been performed at 
low center-of-mass energies ( < 40 GeV) in fixed 
target experiments 0, At y/s ~ 52 — 63 

GeV, the available measurements are not conclu- 
sive due to inconsistencies between different mea- 
surements |l2t 1 1 1| . The measurements at higher en- 
ergy colliders have been at high vr only [15| , or 
have included large uncertainties 0, . Theoreti- 
cal predictions for the RHIC energy region differ sig- 
nificantly (lSl. figj l . Therefore, precise measurements 
of charm cross sections in p+p and d+Au collisions 
in this energy region are crucial. In this paper, 
we report first results on open charm cross sections 
a * -^/s NN = 200 GeV from direct charmed hadron 
D°(D°) reconstruction in d+Au collisions and from 
charm semileptonic decay in both p+p and d+Au 
collisions. These measurements are complementary, 
providing important experimental cross-checks. 

The data used in D° direct reconstruction and 
charm semileptonic decay analysis were taken dur- 
ing the 2003 RHIC run in d+Au and p+p collisions 
a * \/ s nn = 200 GeV with the Solenoidal Tracker at 
RHIC (STAR). A minimum bias d+Au collision trig- 
ger was defined by requiring at least one spectator 
neutron in the outgoing Au beam direction deposit- 
ing energy in a Zero Degree Calorimeter (ZDC). De- 



tailed descriptions of the trigger and centrality def- 
inition in d+Au collisions have been presented in 
a previous publication |2CJ. A total of 15.7 million 
minimum bias triggered d+Au collision events were 
used in the D° analysis. The data samples used in 
the electron analysis in d+Au and p+p collisions 
were described in Ref. |2l| . The integrated luminos- 
ity is about 40 /xb _1 for d+Au collisions and 30 nb _1 
for p+p collisions. 

The primary tracking device of the STAR detec- 
tor is the Time Projection Chamber (TPC) It 
was used to reconstruct the decay of D° — ► K~ir + 
(D° — > K + tt~) which has a branching ratio of 
3.83%. In what follows, we imply (D° +D°) /2 when 
using the term D° unless otherwise specified. The 
exact D° decay topology cannot be resolved due to 
insufficient track projection resolution close to the 
collision vertex. The invariant mass spectrum of 
D° mesons was obtained by pairing each oppositely 
charged kaon and pion candidate in the same event. 
The kaon and pion tracks were identified through 
ionization energy loss (dE/dx) in the TPC wherever 
the identification is possible. Candidate tracks were 
selected having momenta p (px) > 0.3 (0.2) GeV/c 
and pseudorapidity \r/\ < 1. The D° signal with 
Pt < 3 GeV/c and |yj_< 1 after mixed-event back- 
ground subtraction |23| is shown in the left panel 
of Fig. H The signal-to-background ratio (S/N) 
is about 1/600, and the figure of merit (S/y/N) is 
about 6. This distribution was fit to a Gaussian 
plus a linear function to account for the residual 
background not described by the mixed-event spec- 
trum 23] . The open symbols in the left panel of 
Fig. n depict the D° signal after the two-step back- 
ground subtraction. HIJING simulations [24( have 
shown that di-hadron correlations from jets can af- 
fect the line-shape of the background spectrum since 
the shape (slope versus mass) from this contribution 
is different from that of random pairs. To estimate 
the uncertainty in the subtraction of the residual 
background, different normalizations, slopes and fit 
ranges were tried. The resulting uncertainty in the 
D° yield is estimated to be 15%. Within statisti- 
cal uncertainties, the yields of D° and D° are equal. 
The D° — > K~ir + signal could be mis-identified as 
a D° — ► K + tt~ and vice versa when both of its 
daughters are beyond particle identification in the 
TPC. This misidentification results in double count- 
ing which was corrected for in the D° yields through 
a Monte Carlo simulation. 



FIG. 1: (a) Invariant mass distributions of kaon-pion 
pairs from d+Au collisions. The solid circles depict 
the signal after mixed-event background subtraction, 
the open circles after subtraction of the residual back- 
ground using a linear parametrization. (b) dE/dx in 
the TPC vs. particle momentum (p) with a TOF cut of 
|l//3 — 1| < 0.03. Insert: projection on the dE/dx axis 
for particle momenta 1< p <1.5 GeV/c. 



Another detector used in this analysis was a pro- 
totype time-of-flight system (TOFr) jgH based on 
multi-gap resistive plate chamber technology. It cov- 
ers an azimuthal angle A</> ~ 7r/30, and — 1 < r\ < 
0. In addition to its hadron identification capabil- 
ity [23, it allows electrons/positrons to be identi- 
fied at low momentum (j>t < 3 GeV/c) by using a 
combination of velocity information (/3) from TOFr 
and dE/dx measured in the TPC. The right panel of 
Fig-ffldemonstrates th e c i ean separation of electrons 
from hadrons using their dE/dx in the TPC after 
applying a TOFr cut of \l//3 - 1| < 0.03. This cut 
eliminated the hadrons crossing the electron dE/dx 
band. Electrons/positrons were required to origi- 
nate from the collision vertex. Hadron contamina- 
tion was evaluated to be about 10 — 15% in a selec- 
tion optimized for purity and statistics. At higher 
Pt (2 — 4 GeV/c), electrons could be identified di- 
rectly in the TPC since hadrons have lower dE/dx 
due to the relativistic rise of the dE/dx for elec- 
trons. Positrons are more difficult to identify using 
dE/dx alone because of the large background from 
the deuteron band. The hadron contamination in 
this case was found to be < 5% at pr — 2 GeV/c 
and to increase to ~ 30% at px — 3 — 4 GeV/c. This 
was corrected for in the final spectra. Detector ac- 
ceptance and efficiency corrections were determined 
from detailed simulations [2l]]. Total inclusive elec- 
tron spectra from 200 GeV p+p and d+Au collisions 
are shown in Fig. [21 

Gamma conversions 7 — * e + e~ and 7r° — > je + e~ 
Dalitz decays are the dominant photonic sources of 
electron background. To measure the background 
photonic electron spectra, the invariant mass and 
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FIG. 2: Upper panels: Electron distributions from p+p 
(left) and d+Au (right) collisions. Solid and open sym- 
bols depict electrons/positrons ((e + +e~)/2) identified 
via a combination of Time-of-Flight (TOF) and dE/dx, 
and electrons (e~) identified via dE/dx alone. The total 
photonic backgrounds are shown as solid lines. Dashed 
lines depict the various contributing sources. The frac- 
tions were derived from simulations. Bottom panels: the 
ratio of inclusive electrons to the total backgrounds. The 
gray band represents the systematic uncertainty in each 
panel. 



opening angle of the e + e~ pairs were constructed 
from an electron (positron) in TOFr and every other 
positron (electron) candidate reconstructed in the 
TPC 26]. A secondary vertex at the conversion 
point was not required. Simulations with both HI- 
JING |3 and PYTHIA [13 with full detector de- 
scription in GEANT yielded ~ 60% efficiency for 
electrons with pr > 1 GeV/c from such background 
processes. More than 95% of the electrons from 
sources other than heavy-flavor semileptonic decays 
were measured with this method. The remaining 
fraction from decays of 7],uj,p,4> and K was deter- 
mined from simulations. The results are shown as 
solid lines in Fig. The overall uncertainty of the 
background is on the order of 20% and has been in- 
cluded in the systematic errors. Ratios of the inclu- 
sive electrons over the total backgrounds are shown 
in the bottom panels of Fig. [21 The signal is clearly 
in excess of the background above px > 1 GeV/c. 

The non-photonic electron spectra were obtained 
by subtracting the previously described photonic 
background from the inclusive spectra. The re- 
sults are shown in Fig. [21 The D° invariant yields 
d 2 N/ (2-itpTdpTdy) as a function of pr from direct 
reconstruction are shown in Fig. [21 as solid squares. 
Two different fitting methods were used to extract 
dN/dy for the D° at mid-rapidity. In the first 
method, dN/dy was extracted from an exponential 
fit to the D° differential yield in transverse mass 



5 



o 

C\J 





200GeV 


■ 


d+Au (D°+D°)/2 I xZO 






• 


d+Au (e\e)/2 [ x5] 






o 


d+Au e [x5] 






A 

A 


p+p (e + +e')/2 
p+pe' 




"Cr- 








■: z < 







d+Au D° fit [xlO] 








— d+Au c—>e fit [x5j 






:? 


7 /)+/? c— »e fil 








_-— p+pPYTHIA 









12 3 4 



p,- (GeV/c) 

FIG. 3: Reconstructed D° (solid squares) pr distri- 
butions from d+Au collisions at t/snn = 200 GeV. 
Non-photonic electron pr distributions from p+p col- 
lisions (triangles) and d+Au collisions (circles). Solid 
and dashed lines are the fit results from both D° and 
electron spectra in d+Au collisions. The dotted line is 
scaled down by a factor of Ni, in — 7.5 ± 0.4 [20| from 
d+Au to p+p collisions. The dot-dashed line depicts a 
PYTHIA calculation fZft. 



(hit) 23]. In the second method, a simultaneous 
fit was applied to both directly reconstructed Z?°'s 
and the background subtracted non-photonic elec- 
tron distribution in d+Au collisions. For this fit, it 
was assumed that the D° spectrum follows a power 
law in pt from which an electron spectrum was gen- 
erated using the particle composition from |28| and 
the decay generators in PYTHIA. A set of parame- 
ters for the power law was found at the minimum of 
X 2 for the D° and electron spectra. The results are 
shown in Tabled The systematic error is dominated 
by the uncertainties in the background subtraction, 
the extrapolation due to finite pt coverage, and the 
overall normalization (±14% in p+p and ±10% in 
d+Au collisions pcH^I. 

The yield of D° at mid-rapidity is dN/dy = 
0.028 ± 0.004 ± 0.008 and the (p T ) = 1.32 ± 0.08 
GeV/c in d+Au collisions. We used the ratio R — 
N D o /N c5 = 0.54 ± 0.05 from e+e" collider data [H 
to convert the D° yield to a total cc yield. A p+p in- 
elastic scattering cross section of cr^ el = 42 mb was 
used in the calculation and a factor of / = 4.7 + 0.7, 
estimated from simulation 0, |2tJ , was used to con- 
vert the da/dy at mid-rapidity to the total cross 
section. The total charm cross section per nucleon- 
nucleon interaction for d+Au collisions at 200 GeV 
is o»» = dN d D i Au /dy x aZJN^ n Au x f/R = 
1.3 + 0.2 + 0.4 mb from D a alone and 1.4+0.2 + 0.4 
mb from the combined fit of Z?° and electrons. The 
nuclear modification factor |20j was obtained by tak- 
ing the ratio of the electron spectra in d+Au and 
p+p collisions scaled with the underlying nucleon- 
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FIG. 4: Total cc cross section per nucleon-nucleon col- 
lision vs. the collision energy ( ,/8„ „). The dashed line 
depicts a PYTHIA calculation Q. The solid and dot- 
dashed lines depict two NLO pQCD calculations with 
MRST HO, m c = 1.2 GeV/c 2 , fi F = 2m c , fi R = m c and 
2m c , respectively [l^|. 



nucleon binary collisions. It was measured to be 
1.3 ± 0.3 ± 0.3, averaged over 1 < p T < 4 GeV/c. 
This value is consistent with binary scaling within 
the measured errors. 





dN(D°)/dy\ v= o (10- 2 ) 


da? s N /dy\ y =o (mb) 


D° 


2.8+0.4+0.8 


0.29+0.04+0.08 


D°+e ± 


2.9+0.4+0.8 


0.30+0.04+0.09 



TABLE I: dN/dy of D° in d+Au collisions and the corre- 
sponding da/dy of cc pair per nucleon-nucleon collision 
at ,/sZZ= 200 GeV. 



The beam energy dependence of the cross section 
is shown in Fig. H Both default PYTHIA [27| and 
NLO pQCD [18j calculations reasonably describe the 
results at lower energies, but underpredict the total 
charm cross section at y/s NN — 200 GeV. A NLO 
pQCD calculation (solid line) with fragmentation 
and renormalization scales chosen to be /if = "2.Tn c 
and /xr = m c (m c = 1.2 GeV/c 2 ) reproduces our re- 
sult. The underprediction by PYTHIA of the charm 
cross section is also evident in Fig. [3] , the charm de- 
cayed electron px distribution shown as dot-dashed 
line. Furthermore, the slope of the PYTHIA distri- 
butions is much steeper than the measured distribu- 
tion. There are also indications that a large charm 
production cross section at ^s NN ~ 300 GeV is es- 
sential to explain available cosmic ray data |2flj |. 

At RHIC energies, binary scaling of the open 
charm production is expected between p+p, p+A 
and A+A collisions Q. If correct, the results of 
this study suggest a much larger charm yield in cen- 
tral Au+Au collisions than previously assumed in 
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statistical thermal models @, IS @ based on some 
pQCD/PYTHIA calculations. This would rule out 
several predictions [3, 13, |9j of charm production not 
previously excluded by the upper limit (below binary 
scaling) set by J/ip production in central Au+Au 
collisions [3(j ■ Future heavy ion runs at RHIC with 
open charm and J ftp measurements will enable us to 
study the flow and thcrmalization of charmed parti- 
cles. 

In summary, the charm cross section and trans- 
verse momentum distribution for p+p and d+Au 
collisions at ^/s NN = 200 GeV have been measured 
by the STAR collaboration at RHIC. Independent 
measurements of the reconstructed D° and single 
electrons from charm semileptonic decay are con- 
sistent. The total cross section at this energy was 
compared to theoretical calculations. The result has 



important consequences for charm quark coalescence 
in Au+Au collisions at RHIC. 
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